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In this paper we present an electromagnetic modal theory for 
characterizing parabolic-index multimode fiber splices with either 
intrinsic or extrinsic mismatches. The theory agrees with previously 
published theoretical results for transverse offset using a uniform 
power distribution. It also agrees with new experimental measure- 
ments made with a long, spliced input fiber using a published, 
theoretical, steady-state modal power distribution. This modal theory 
predicts, and experiment confirms, a previously unreported periodic 
fluctuation in splice loss as a function of wavelength for intrinsic 
parameter mismatch. The analysis also predicts a large degree of 
mode mixing for transverse offset but negligible mode mixing for 
parameter mismatch in typical multimode fiber splices. 

I. INTRODUCTION 

Theoretical predictions of splice loss in multimode fibers have been 
attempted by many researchers (see Ref. 1 for a list). However, these 
analyses 1-3 do not adequately predict measured splice loss results, and 
most do not address mode mixing effects, which are important for the 
prediction of the bandwidth of concatenated lengths of fiber. The best 
agreement with splice loss measurement data is an empirical model 
based on geometric optics. 2 However, such an analysis of splice loss 
provides only a limited description of the effect of splices. Electromag- 
netic theory, using the coupling coefficients of individual modes, pro- 
vides a complete treatment of splice loss and mode mixing. The only 
published electromagnetic theory 3 appears to be in error, since predic- 
tions do not approach the well-known correct geometric optics limit 
for splice loss of fibers with a large normalized frequency, V, and a 
uniform modal power distribution. This paper presents an electromag- 
netic analysis for splices, which gives single-term expressions for mode 
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coupling coefficients for the case of either transverse offset or profile 
parameter mismatch for parabolic-index multimode fibers. These cou- 
pling coefficients are used to calculate loss and degree of mode mixing 
in splices. The results approach the proper geometric optics limit for 
a uniform modal power distribution. Furthermore, if we compare 
theoretical and experimental loss to transverse offset results, we can 
verify a published theoretical "steady-state" power distribution 4 in 
multimode parabolic-index fibers. Comparing the results of this theory 
with a measurement of splice loss as a function of transverse offset 
gives the modal power distribution in a fiber. A wavelength dependence 
of splice loss for intrinsic parameter mismatch is predicted by this 
theory and has been experimentally verified. Because splices change 
the modal power distribution, the power redistribution can cause 
additional loss as the system evolves towards the steady state in the 
receiving fiber. In the past 2 this loss has been considered part of the 
total splice loss, but the improvement of fiber quality has reduced this 
effect significantly for typical lengths between splices. Because we 
have ignored the power redistribution effects (after the splice), these 
theoretical results are only strictly valid for short lengths of fiber after 
the splice; however, they should remain valid for typical distances 
between splices. 

II. THEORY 

The modal amplitude coupling coefficient, C\m\a.\\ 7112(12], for mode 
m.\a\ of the transmitting fiber and mode m.2ci2 of the receiving fiber is 
obtained directly from the theory of excitation of weakly guiding 
fibers. 5 (The variable m is the radial mode number and a is the 
azimuthal mode number.) 

C[miai;7w 2 a 2 ] =-( — J \ Em^E^dA, (1) 

where for parabolic-index fibers (using the infinite profile approxima- 
tion) the electromagnetic modal field, E miBll , propagating in the z 
direction in the fiber core is: 

*** = (f) A„J-TLW(t)e~* r^'*, (2) 

where x and y are the linearly polarized unit field vectors of the mode. 
e co and /* are the permittivity and permeability of the core, respectively, 
and € co = ricoe. Then 
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(3) 



and 



a 



^ 2 = (x"")- 2iV ^ (4) 

V= — n co aV2A, 

A 

where the principal mode number is N(= 2m + a + 1), and Lm(t) is an 
associated Laguerre polynomial. The variable a is the fiber core radius 
and A = (rico — n c \)/n co is the normalized maximum refractive index 
difference. For weakly guiding fibers the polarization property of the 
incoming mode is preserved by the splice, as can be seen from 

eq. (1). 

In general we can write: 

C(miai; m 2 a 2 ) = K-I(mion; m 2 a 2 ), (5) 



where 



K = 



m\\m 2 \ 



Tl/2 



(mi + |ai|)!(m2+|a 2 |)!_ 
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7(miai ; m 2 ct 2 ) is a function that depends on the type of splice mismatch. 
Table I shows the expressions for I(micti,m 2 a 2 ) for the two cases of 
intrinsic parameter (V) mismatch and transverse offset between iden- 
tical fibers (r„). The derivation of these equations is given in the 
appendix. 

For the parameter mismatch case, the coupling coefficients are a 
function of the ratio of the normalized frequencies of the two fibers. 
The hypergeometric function, 2 Fi, is, in this case, a power series in y 2 
of order mi . Note that because azimuthal symmetry is preserved in 
the splice, only modes with the same azimuthal mode number couple. 
In the identical fiber transverse offset case, the coupling coefficients 
are simple products of Laguerre-Gaussian polynomials with argument 
proportional to the normalized offset (r /a). In this case all azimuthal 
modes have finite coupling coefficients, but, for small offsets, coupling 
is much stronger for nearest neighbor azimuthal modes. 

Assuming a random phase relationship between the modes of the 
transmitting fiber, the total power coupled into mode m 2 a 2 , Pm 2 <, 2 , is: 



n^c 



T 

= X 



\C(miai\ m 2 a2)| P> 



m, ai , 



(7) 



and the total splice loss, S s , is then: 
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The mode mixing factor, T (i.e., the fraction of the power that is 
redistributed into different modes), is: 

T 

£ |C(mia 1 ;m 2 « 2 )| 2 P miai 



r = i- 



" , l =,n 2' a l =a 2 



R 



(9) 



The superscripts R and T refer to sums over the bound mode spectra 
of the receiving and transmitting fibers, respectively. 

The range of allowed radial and azimuthal mode numbers, m and 
a, for a fiber with normalized frequency V is: 



m *S m max - INT(V/4) 
a<<w = INT(V/2)-l 

Ntns* = INT(V/2), 

Table I — Mode coupling functions 
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where INT indicates the integer part. It is essential to realize the 
discrete nature of the bound mode spectra to understand the splice 
loss predictions for intrinsic parameter (V) mismatch. The explicit 
wavelength dependence of splice loss and mode mixing is demonstrated 
by the presence of V in the coupling coefficients. However, an implicit 
wavelength dependence due to the allowed spectrum of bound modes 
occurs from eq. (10). 

III. RESULTS 

Equations (8) and (9) can be used to calculate splice losses and 
degree of mode mixing for the two cases of parameter mismatch and 
transverse offset. To validate this analysis, they must be compared 
with existing theories wherever possible. Most existing results are 
geometric optics evaluations of splice loss, and the modal power 
distribution equivalent to any assumed ray distribution is difficult to 
evaluate in general. However, the ray distribution equivalent to the 
uniform modal power distribution has been well documented, 6 and 
comparison of splice loss predictions for this case demonstrates that 
this analysis is correct. Calculations and experimental results for the 
effects of different modal power distributions on the mode mixing and 
loss at a splice are shown in the sections that follow. 

3.1 Identical fiber transverse offset 

This section presents the results for mode mixing and splice loss 
associated with transverse offset of the fiber axes for identical fibers. 

3.1.1 Splice loss 

Splice loss predictions for transverse offset with a uniform power 
distribution for several values of the normalized frequency, V, are 
shown in Fig. 1. It is well known that, as V increases, the electromag- 
netic analysis should asymptotically approach the geometric optics 
predictions. This is confirmed in Fig. 1. The only other published 
electromagnetic splice loss theory 3 does not appear to converge to the 
correct geometric optics limit, e.g., from Fig. 3 of Ref. 3, the splice loss 
for 0.2a offset (for 2a = 50 /an) is -1.0 dB, whereas geometric optics 
predict a value of 0.8 dB. The small deviation of our results for small 
offset (<0.1a) is believed to be caused by coupling of power to the 
bound modes of the receiving fiber through the evanescent cladding 
fields. This effect cannot be predicted using ray optics. 

Uniform- power distribution results do not agree with realistic splice 
loss measurements. 2 Figure 2 demonstrates loss versus transverse 
offset for several power distributions, the Gaussian splice loss model, 2 
and experimental data. 8 The experimental data shown were obtained 
using 0.82-/nm laser excitation of a 7-km input fiber containing about 
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Fig. 1 — Splice loss vs. transverse offset with uniform input power. 

14 splices into a final splice that was offset in both x and y orthogonal 
axes. The offset data are in good agreement with previous data and 
should represent a realistic "steady-state" condition. The theoretical 
power distribution for steady state due to microbending 4 gives excellent 
agreement with this data. (A power distribution that has been used to 
approximate this steady-state distribution is also plotted to show the 
sensitivity of the loss to the choice of power distribution.) The excellent 
agreement between the data and this analysis using the steady-state 
power distribution confirms both the distribution and the splice loss 
theory. Although this power distribution gives excellent agreement 
with the data for very long lengths of fiber before the splice, it does 
not agree as well as the Gaussian model 2 with shorter-length input 
fiber measurements because the power distribution has not achieved 
"steady state." 

Figure 3 compares experimental data for shorter lengths of input 
fiber (~1 km) and different excitation conditions 8 with theoretical 
predictions using both the Gaussian model 2 and various power distri- 
butions in this modal analysis. Note that the different excitation 
conditions significantly affect the splice loss for larger offsets. The 
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Fig. 2— Loss vs. transverse offset for different modal power distributions. 

empirical model of Ref. 2 agrees well with loss measurements for 
typical offsets expected in practice of 0.1 to 0.2 (r /a) for both sources 
and agrees very well over a wider range for the laser source. Modal 
power distributions for the electromagnetic analysis to describe this 
situation were empirically chosen so that, for very long fiber lengths, 
they would degenerate to the steady-state distribution. The simplest 
choice for a power distribution satisfying this criteria is: 



P(miai) = Jo 



2.405 



N 



N n 



(1 - e" L ) 



(11) 



The choices of tj shown in Fig. 3 (for L = 1 km) demonstrate that 
excellent agreement between the modal analysis and the splice loss 
versus transverse offset for the shorter input fiber lengths is possible 
for both sources. Furthermore, for rj = 3, the results agree with the 
predictions of the empirical model 2 within 0.03 dB. The theoretically 
predicted steady-state distribution gives optimistic results for this 
short length. The existence of a modal power distribution that is in 
good agreement with the results of Ref. 2 for transverse offset supports 
the basic assumptions of the empirical model. 

3.1.2 Mode mixing 

Equation (9) can be used to calculate the amount of power in the 



FIBER SPLICES 1195 



1.2 



0.6 - 



0.4 - 



A 


0.82 jllm LASER SOURCE 


O 


0.82 )xm LED SOURCE 
- GAUSSIAN MODEL (REF. 2) T 





- STEADY STATE - J [2.A{X\ V2 ] (REF. 4) > 

- J [2.4 (l-e^X'*] / 








-7 [2.4(1-e- 25 )X 1/2 ] / 




°/A 




/'/ 




/'/ / 


- 


/// 




// / 




// / 




V / 




/// 






// / 








'yy/ 


'Cf 


jj>^ 


^2^5 


i i i 



0.1 0.2 0.3 0.4 

NORMALIZED TRANSVERSE OFFSET IN r a la 

Fig. 3 — Splice loss vs. transverse offset for 1-km input fiber. 



receiving fiber that has changed propagation mode because of the 
splice. By changing the summation criteria one can calculate any 
desired mode mixing factor, such as power coupled within the same 
principal mode group, nearest neighbor mode groups, etc., where the 
summation is over all mi, ol\ and m.2, a.2 that differ in principal mode 
number by AiV, which is written as: 



2 for A2V = 0, 1, 2, 

2m2+a 2 -2m,-a,=A7V 



(12) 



Mode mixing as a function of normalized transverse offset for two 
different input power distributions is shown in Fig. 4 as the percentage 
of receiving fiber power that has been redistributed into neighboring 
mode groups. These differ in principal mode number, N, by AN of 0, 
1, 2, 3, and 4. (AiV = represents power coupled within the same 
degenerate mode group in both transmitting and receiving fibers.) The 
power distributions used were the uniform and the theoretically pre- 
dicted power distribution for the steady state for microbending. 4 It is 
important to note that the percentage of the power-changing propa- 
gating mode groups is very high even for small offsets with low splice 
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Fig. 4 — Power coupling into neighboring mode groups for transverse offset. 

loss, i.e., 31 percent for 0.04a offset and 0.07 dB splice loss for a V = 45 
fiber with the steady-state power distribution. Mode mixing is a fairly 
weak function of the modal power distribution, which indicates that 
the mode mixing is uniform over all the mode groups as can be seen 
from the local numerical aperture arguments. 2 As offset increases, the 
mode mixing changes from being primarily to the nearest neighbor 
mode groups (AN = ±1) to being redistributed over a wider range of 
principal mode groups. The strength of the mode coupling for even 
small offsets with low splice loss, and the relatively small offsets at 
which it spreads over many mode groups, are both initially surprising. 

3.2 Parameter mismatch 

For fibers with no axial offset, the individual mode fields are insen- 
sitive to small changes in the V parameter, so that for splices with 
only slight parameter mismatch, as can be seen in eq. (1), self-coupling 
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dominates (i.e., mode m,i<x\ of the transmitting fiber couples primarily 
with mode miai of the receiving fiber). For example, there is less than 
3-percent mode mixing for a 5-percent parameter mismatch with 
uniform modal excitation over the wavelength range of 0.8 to 1.4 /mi. 
Therefore, splice loss due to parameter mismatch is caused only by 
this slight field mismatch unless the ranges of the bound mode spectra 
of the two fibers [eq. (10)] differ, i.e., if the normalized frequencies of 
the fibers are such that ANmax > 1. In that case: 

AN™, = NL*-NL* (13) 



or 



AAU* - INT 



y-n co {a T ^2k~T) 



-INT 



t- n co (a r V2Ar) 



(14) 



When AiV ma]t is zero, the splice loss is small, as discussed above. 
However, when AiVma* > 0, the modes in the highest-order mode group 
of the transmitting fiber are not bound modes in the receiving fiber. 
Because self-coupling dominates, most of the power in the highest- 
order modes couples to leaky /radiation modes in the receiving fiber 
and is therefore lost. The large degeneracy of the higher-order mode 
groups [eq. (4)] accentuates this loss. From this discussion we see that 
parameter mismatch splice loss is caused primarily by the difference 
between the bound mode volumes of the two fibers, which is a function 
of the normalized frequencies, V. Therefore, both A and radius mis- 
matches can be expressed simply as V mismatch. (An x-percent radius 
mismatch is equivalent to a 2*-percent A mismatch for small degrees 
of mismatch.) 

Theoretical predictions of splice loss versus normalized parameter 
mismatch (AV/V) for uniform and steay-state input power distribu- 
tions at 0.82 and 1.3 jam are shown in Fig. 5. The discontinuities of the 
splice loss curves are caused by the discreteness of the bound mode 
spectra, as can be seen from the definition of AiVmax in eq. (14). In Fig. 
5a this discreteness is most obvious due to the choice of a uniform 
input power distribution, which accentuates the power in the highest- 
order mode group. Note that at the longer wavelength the loss is 
initially higher for a given parameter mismatch because the percentage 
of the total power that is in the highest-order mode groups is greater 
than at shorter wavelengths. Each discontinuous increase of splice loss 
is caused by the stripping of an additional high-order mode group of 
the input fiber, which does not propagate in the receiving fiber as 
A V/V increases. The splice loss predicted by geometric optics (which 
ignores the discreteness of the mode spectra) is also shown. 1 The 
results for the steady-state power distribution used previously are 
shown in Fig. 5b displaying similar discontinuous features. However, 
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Fig. 5 — Splice loss vs. normalized parameter mismatch in &V/V for (a) uniform power 
distribution and (b) steady-state power distribution. 

the decrease in the order of magnitude of the loss allows the observa- 
tion of the small loss component caused by the slight field mismatch 
between similar modes of the two fibers. This is reflected in the small 
slope of the levels. Because the steady-state distribution assumes that 
there is no power in the highest-order mode group of the transmitting 
fiber, the first discontinuous change of splice loss of Fig. 5a is negligible 
in Fig. 5b. 

Figure 6 shows the percentage of mode mixing at a splice caused by 
parameter mismatch for 0.82- and 1.3-jum wavelengths, and the uniform 
and steady-state power distributions. The substantially higher mode 
mixing for the uniform power case compared to the steady-state 
indicates that mode mixing is dominated by the higher-order modes. 
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Fig. 6 — Percent of mode mixing vs. normalized parameter mismatch. 

This contrasts with the transverse offset ease shown in Fig. 4, where 
mode mixing is essentially the same for uniform and steady-state 
distributions. Reduced mode mixing at longer wavelengths is due to 
the more diffuse bound mode fields, which result in stronger self 
coupling between modes. In comparison to the transverse offset case, 
the degree of mode mixing for the case of parameter mismatch is 
substantially reduced. Furthermore, a detailed analysis of the mode 
mixing shows that the coupling is almost completely nearest neighbor 
coupling (AN = 1), e.g., for AV/V ~ 0.2, 38 percent of the total 44- 
percent mode mixed power is coupled to nearest neighbor modes. The 
deviations from a smoothly increasing function in these curves occur 
at the transition points of the splice loss curves (see Fig. 5). This is 
caused by the transition of the highest-order mode group in the 
receiving fiber from propagating to lossy, removing its contribution to 
the overall mode mixing. 

3.3 Wavelength dependence of splice loss for parameter mismatch 

As we can see from eq. (13), for a given degree of parameter 
mismatch, ANmax is also a discontinuous function of wavelength. For 
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example, if, at some particular wavelength, Vt is 46.1 and Vr is 45.9, 
ANmax is 1 and splice loss is -0.5 dB. A slight change of wavelength 
(~1 percent) will cause both V's to simultaneously be greater or less 
than 46 with AiV,^ = 0, and splice loss is -0.02 dB. Theoretical 
considerations [eq. (14)] also indicate that the period of the splice loss 
fluctuations should increase with increasing wavelength because the 
mode volume is proportional to 1/A. The magnitude of the fluctuation 
should also increase with wavelength because the relative amount of 
power in the highest-order mode group increases as the mode volume 
decreases. The splice loss wavelength dependence for a parameter 
mismatch (AV/V) of 0.04 and uniform input power distribution is 
shown in Fig. 7. The magnitude of the splice loss changes considerably 
over the entire wavelength range. Splice loss calculated from geometric 
optics for this parameter mismatch case is also shown in Fig. 7 and is 
approximately the wavelength averaged splice loss as calculated from 
modal theory. Although the splice loss displays this pathological 
behavior, the mode mixing remains relatively small and continuous. 

Figure 8 gives two further examples of the theoretical splice loss 
versus wavelength for uniform power excitation and core radius mis- 
matches of 1.0 and 1.5 /mi. For Aa of 1.5 jum, there are ranges of 
wavelength over which AiVmax — 2, resulting in even larger fluctuations. 
The choice of a uniform power distribution unrealisticalry enhances 
the magnitude of this effect. For example, with the choice of the 
steady-state distribution, there would be very little power in the 
highest-order mode group, and therefore this change in loss would be 
very small, as we can see in Fig. 5b. Splice loss is not expected to vary 
as abruptly as shown in Figs. 7 and 8, even for a uniform power 
distribution, because the cladding splits the degeneracy of the highest- 
order mode group and some modes of this group remain bound even 
when the group itself, as predicted by the infinite parabolic approxi- 
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Fig. 7 — Splice loss and mode mixing vs. wavelength for A V/V = 0.04. 
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mation, is not bound. This results in a smoothing of the loss versus A 
curve, but the periodicity should not be affected. The observation of 
this effect would further verify this analysis. 

3.4 Experiment 

The experimental setup shown in Fig. 9 was designed to enhance 
the power in higher-order mode groups in an attempt to observe the 
predicted wavelength dependence of splice loss shown in Figs. 7 and 8. 
A 5 m, 63-ju.m core, step-index fiber was over-filled and used to excite 
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Fig. 9 — Splice loss measurement system. 
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Fig. 10 — Measured period and splice loss for parameter mismatch. 



the test fiber, which was wound on a 7.5-cm mandrel. A cut-back test 
determined that approximately 15 m of fiber on this mandrel success- 
fully eliminated leaky modes in the transmitting fiber. This fiber was 
spliced with micropositioners to a 1 to 2 m length of fiber, which 
differed in core radius by 1.0 ± 0.2 /xm (but with the same A), as 
determined by refracted near-field measurements. Monochromator 
resolution was 2 nm and the repeatability of the splice loss measure- 
ments was 0.02 dB. As a final check the test fiber was spliced to a 
length of identical fiber. Maximum splice loss change seen in this case 
was 0.02 dB, with no discernible wavelength dependence. 

The measured wavelength dependence of splice loss and its period 
are plotted in Fig. 10 for the mismatched fibers. The data show 
excellent agreement with the period versus wavelength curve obtained 
from Fig. 8. The increase in maximum amplitude with wavelength seen 
in Fig. 10 is in qualitative agreement with the theory. Cladding effects 
(ignored in the infinite parabolic approximation) and only partial 
success in achieving a uniform power distribution are the probable 
reasons that the maximum amplitude is less than expected from 
Fig. 8. 

IV. CONCLUSIONS 

In this paper we have presented a modal theory that allows the 
calculation of mode mixing and loss resulting from parabolic-index 
multimode fiber splices. The predictions of this model have been 
confirmed by experiment, including a previously unreported wave- 
length dependence of splice loss for intrinsic parameter mismatch. 
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This phenomenon is not expected to be significant for practical 
"steady-state" system losses but is important for providing insight into 
the actual loss mechanisms at a splice as well as verifying the validity 
of the theory. Another effect of this discreteness of the bound mode 
spectrum has previously been reported for the wavelength dependence 
of bending loss. 9 The theory also successfully describes splice loss 
versus transverse offset phenomena and has been used to verify the 
predicted steady-state power distribution. The analysis provides a 
basis for understanding mode-mixing effects at splices, which is im- 
portant for bandwidth studies on concatenated fiber lengths. Splice 
losses caused by offset or parameter mismatch produce greatly differ- 
ent degrees of intermodal coupling. This implies that the effect of the 
splice on the bandwidth of a concatenated length is insignificant for 
the case of small parameter mismatch, but may be significant for 
transverse offset (due to the strong mode mixing). The sensitivity of 
the splice loss predictions to the choice of the input power distribution 
may allow this analysis to provide an evaluation of modal power 
distributions in fibers from experimental measurements of splice loss 
versus transverse offset. 

The particular case of longitudinal offset has been obtained from a 
generalization of the parameter mismatch theory but has not been 
treated here due to its relative unimportance in practical systems. 
Gloge used geometric optics to show that fiber axis tilt at the splice is 
equivalent to a transverse offset via the relationship (r /a) 2 is equal to 
sin 2 0/2A, where 6 is the angle between the fiber axes. 10 This equivalence 
is also valid from the electromagnetic analysis for pure tilt or offset. 
The modal analysis of combinations of offset (both transverse and 
longitudinal) with tilt and parameter mismatch (as discussed in Ref. 
10) are subjects of continuing work. 
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APPENDIX 

Derivation of Modal Coupling Coefficients 

The principal steps in the derivation of the modal coupling coeffi- 
cients (eq. 1) presented in Table I are shown in this appendix. 

A.1 Intrinsic parameter (V) mismatch 

In this case both fibers have fields as defined in eq. (2), where the 
argument of the field, Z, is 

Z T =V T (r/a T ) 2 (15) 

and 

Z R - Z T (1 + e), (16) 

where 

V T = (27r/X).n C o-y/2K^-a T (17) 

and 

e-^M|-l. (18) 



al \V 



T 



where the normalization expressions A ma also contain Vt and Vr . 

Examination of the azimuthal integral of eq. (1) shows that, because 
there is no transverse offset of the fiber axes, azimuthal symmetry is 
preserved and only modes with the same azimuthal symmetry couple. 

Changing the radial integral from an integration over r to one over 
Zt, we find that it is in the form of a well-documented integral of 
Laguerre Gaussian functions, 11 giving the expression shown in Table I. 
The hypergeometric function, 2 Fi(—m 2 , -mi, -(mi + m 2 + a), y 2 ) is a 
finite power series in v 2 whose order is the lesser of mi or m 2 , since 
both mi and m 2 are integers. 
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A.2 Transverse offset— identical fibers 

The modal fields to be used in eq. (1) can still be expressed by eq. 
(2) for both fibers. However, the radial and azimuthal variations of the 
receiving fiber {tr, 4>r) need to be expressed in terms of the transmit- 
ting fiber (rr, <>) and the transverse offset, r Q , of the axes. Using 

r R = rr + r 2 - 2r T r cos </> (19) 

and the integral relation for Laguerre-Gaussian functions 12 

(—1) m f" 

e -z/2 Z la|/2 L la| (Z)= !_L j a (^xZ)e-^ 2 x^ 2 L^(x)dx (20) 

Jo 

we can then use the addition theorem for circular cylindrical waves 13 
to write the receiving fiber field amplitude £«(m, a) as 

E R (m,a) = { -^f- i e- i{p+a) * 



p — a 



f 



where 



J P (y/xYo)J p+a (y/xZ^)e- x/2 x^ /2 L^(x)dx, (21) 

Z = V(r /a) 2 

Z T = V(r T /a) 2 . (22) 

The azimuthal integral of eq. (1), for coupling between mode (mion) 
of the transmitting fiber and mode (7*120:2) of the receiving fiber, reduces 
the infinite summation over/), introduced by eq. (21), to a single term 
for 

p + a2 = ai. (23) 

After a change of variable from rr to Zt defined above, the radial 
integral of eq. (1) can be evaluated by again using eq. (20). (In 
evaluating this integral it is imperative to realize that the radial field 
variation of the modes is a function only of the magnitude of the 
azimuthal mode number, \a\.) This evaluation produces two cases, 
which must be treated separately: 
(i) cti and 02 — both the same sign 

(ii) oti and 0:2 — different sign. 

For both cases the final integral over the dummy variable x, intro- 
duced by the Laguerre-Gaussian integral identity [eq. (20)], can be 
evaluated by using the relationship 14 

Ln s (x) = (-l) s x s ^—^LLs(x) (24) 
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to write the integral in the form 

/= f y^e- ay2 L p m -°(ay 2 )L° n (oty 2 )J P (x-y)dy, (25) 



which is a well-documented integral giving 

/ = (-l) m+n (2a)- (p+1) x p e- x2/4a L° m - m+n (x 2 /4<x) .L p -° +m - n (x 2 /4a). (26) 

In both Refs. 15 and 16, the formula is slightly incorrect due to a 
reversal of the subscripts of the Laguerre polynomials. The formula of 
Refs. 15 and 16 is seen to be incorrect by noting that eq. (25) above is 
a Hankel transform of order p, so that, by applying the inverse 
transform to the result, one should expect to obtain the original 
function, i.e., the argument of the original Hankel transform. The error 
persists in Ref. 15 because it used Ref. 16 as its source for this result. 
The use of this relationship then provides the results in Table I. 
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